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Abstract. The purpose of this preliminary work is to determine if Large Magellanic Cloud 
(LMC) Cepheids have stellar winds. If a Cepheid undergoes mass loss then at some distance 
from the star, a fraction of the gas becomes dust, which causes an infrared excess. Mass loss 
is tested using OGLE II optical observations and SAGE infrared observations for a sample of 
488 Cepheids. The resultant mass-loss rates range from 10" 12 to 10" 7 M e /yr. Using the mass- 
loss model we compute infrared stellar luminosities for the sample of Cepheids and compare 
predicted infrared PL relations with observed relations. The predicted relations not only vary 
from the observed relations, implying mass loss plays a significant role, but also show evidence 
for non-linearity. It is determined that mass loss is important for LMC Cepheids. 
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1. Introduction 

LMC Cepheids are important tools for testing stellar astrophysics and for calibrating 
the Period-Luminosity (PL) relation at various wavelengths. Recently infrared PL rela- 
tions have been determined for LMC Cepheids (Ngeow & Kanbur 2008; Freedman et al. 
2008). These PL relations are more universal because the pulsation amplitude decreases 
at longer wavelengths and metallicity contributes less. 

Infrared observations of Galactic Cepheids have found the existence of circumstellar 
shells (Kervella et al. 2006; Merand et al. 2006; Merand et al. 2007) and it has been 
suggested that the shells are caused by stellar winds with dust shells forming at large 
distance from the star (Neilson & Lester 2008). The dust causes the observed infrared 
excess. Mass loss is also a potential solution to the Cepheid mass discrepency which 
is the difference between mass estimates from stellar evolution calculations and stellar 
pulsation models (Cox 1980). If LMC Cepheids have significant stellar winds then mass 
loss can explain the discrepancy and may also cause infrared excess. 

2. Method 

We wish to determine mass-loss rates for LMC Cepheids. To do this we use BV and I 
observations from the OGLE II survey and IRAC observations from the SAGE survey. 
Having up to seven fluxes as a function of wavelength, we can x 2 fit mass-loss rates and 
radii by predicting stellar luminosities and circumstellar shell luminosities, which are 
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Figure 1. The apparent magnitudes LMC Cepheids at infrared wavelengths. 
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The effective temperature of a Cepheid is determined using the Temperature-Period- 
Color relation from Beaulieu et al. (2001). The dust grains are assumed to have size a 
that ranges from 0.005 to 0.025 /zm, and the density of a grain is about lg/cm 3 , and 
is the absorption efficiency. The mass-loss rate is that of the dust and dust velocity 
is about the escape velocity of a Cepheid, 100km/ s. The quantity W(r) is the dilution 
factor and is a function of stellar radius. The dust temperature scales as the effective 
temperature. Therefore given a dust mass-loss rate and radius we can predict the total 
luminosity (stellar + shell) of a Cepheid at various frequencies and compare these to 
observations. The gas mass-loss rate is assumed to be 250 times larger than the dust 
mass-loss rate but this choice of dust-to-gas ratio leads to a minimum prediction of gas 
mass-loss rate. We apply this analysis to a sample of 488 Cepheids with observed infrared 
fluxes shown in Figure 1 (Left). 



3. Results 

The gas mass-loss rates and best-fit \ 2 values are shown in Figure 2. The values 
imply that mass loss is important for LMC Cepheids and contribute to the observed 
infrared fluxes. We take the predicted stellar fluxes and compute infrared PL relations 
to understand the contribution of mass loss on the observed relations. The relations are 
shown in Figure 3 (Left) and the slopes and zero points are given in Table 1. 

The infrared PL relations predicted here differ from those found in the literature 
(Ngeow & Kanbur 2008; Freedman et al. 2008) as can be seen in Figure 3 (Right). The 
PL relations from Ngeow & Kanbur (2008) tend to have brighter zero points, and sim- 
ilar slopes, resulting from larger infrared excesses from the shortest period Cepheids. 
The predicted stellar IR PL relations have shallower slopes than the relations from 
Freedman et al. (2008)|but similar zero points. 



It has been shown that LMC PL relations are non-linear in the optical wavelengths 
(Ngeow et al. 2005) with the K-band relation being marginally linear (Ngeow & Kanbur 
2006). F-tests of the Ngeow fc Kanbur (2008)| found that the IR PL relations were linear 
from 3.6 to 5.8 fim while the 8.0 fim relation is non-linear. The optical non-linear PL 
relations are to have a period break at 10 days, where the slope is steeper for P < lOd 
than for P > lOd. The 8.0 /im relations slope is more shallow for the P < lOd. We 
apply the F-test (Kanbur & Ngeow 2004) to the predicted stellar fluxes and find that the 
predicted IR PL relations are consistent with being non-linear in a similar fashion as the 
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Table 1. Best Fit Parameters for Predicted PL Relations 



Type 




Slope 


Zero Point 


Dispersion 


Linear 


3.6 


-3.14 ± 0.024 


15.99 ± 0.017 


0.110 




4.5 


-3.15 ± 0.023 


15.91 ± 0.017 


0.108 




5.8 


-3.16 ± 0.023 


15.92 ± 0.017 


0.107 




8.0 


-3.17 ±0.022 


15.92 ±0.016 


0.105 


Non- 


3.6 


-3.24 ± 0.038 


16.04 ± 0.025 


0.107 


Linear 


4.5 


-3.24 ± 0.037 


15.97 ± 0.024 


0.106 


P < Wd 


5.8 


-3.25 ±0.037 


15.97 ±0.024 


0.104 




8.0 


-3.26 ±0.036 


15.97 ±0.024 


0.103 


Non- 


3.6 


-2.97 ±0.13 


15.81 ±0.15 


0.125 


Linear 


4.5 


-2.99 ±0.12 


15.75 ±0.14 


0.122 


P > Wd 


5.8 


-3.00 ±0.12 


15.75 ±0.14 


0.120 




8.0 


-3.02 ±0.12 


15.76 ±0.14 


0.118 




0.4 0.6 0.8 1 1.2 1.4 1.6 0.4 0.6 0.8 1 1.2 1.4 1.6 

Log P(d) Log P(d) 

Figure 2. (Left) The values of \ 2 f° r the fits to the sample of Cepheids. (Right) The predicted 
mass-loss rates of the Cepheids. The circled points represent Cepheids where the infrared ob- 
servations are likely blended stars. 

optical PL relations wit the 8.0 /im data being marginally linear. The slopes and zero 
points of the non- linear relations are given in Table 1 and the comparison of the predicted 
linear and predicted non-linear relations are shown in Figure 3 (Right). It appears mass 
loss masks the non-linearity by causing a larger infrared excess in the shorter period 
Cepheids than in the long period Cepheids. This is because the long period Cepheids are 
already redder and the mass-loss rates on average are similar for all periods. This also 
explains the non- linear behavior of the observed 8.0 /im relation and explains why the 
K-band relation is marginally linear. 



4. Conclusions 

In this work, we have modeled mass loss in Cepheids as a spherically symmetric wind 
where dust grains form at some distance from the star where the temperature of the 
wind decreases to the condensation temperature. The dust in the wind causes an infrared 
excess that can be observed. By applying this model to OGLE BVI and SAGE IRAC 
observations we can determine mass-loss rates and analyze the role mass loss plays. We 
find the following conclusions: 

• Mass loss is significant in LMC Cepheids ranging from 10~ 12 to 10~ 7 M Q /yr. 

• The mass- loss rates imply that Cepheids may lose up to 1M on the second cross- 
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Figure 3. (Left) The predicted apparent magnitudes of the sample of LMC Cepheids at infrared 
wavelengths. (Right) The comparison of the observed IR PL relations to our predicted linear 
PL relation. 

ing of the instability strip which is consistent with the mass discrepancy (Keller 2008). 
Therefore mass loss is a plausible solution. 

• Infrared excess due to mass loss have a significant contribution to the IR PL relations 
affecting the slope and zero point. 

• The predicted stellar IR PL relations show evidence for non-linearity. 

• The computed mass-loss rates are dependent on the dust-to-gas ratio which is metal- 
licity dependent. Because mass loss affects the structure of IR PL relation then the 
observed IR PL relations are metallicity dependent. 
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